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l/-) ■ We present new investigation of the Lamb shift (2P]/ 2 — 2S 1 / 2 ) in muonic deu- 

terium {lid) atom using the three-dimensional quasipotential method in quantum 
electrodynamics. The vacuum polarization, nuclear structure and recoil effects are 
calculated with the account of contributions of orders a 3 , a 4 , a 5 and a 6 . The results 
are compared with earlier performed calculations. The obtained numerical value 
of the Lamb shift 202.4139 meV can be considered as a reliable estimate for the 

J; comparison with forthcoming experimental data. 
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The muonic deuterium (/j,d) is the bound state of negative muon and deuteron. The 
lifetime of this simple atom is determined by the muon decay in a time r M = 2.19703(4) • 10 -6 
s. When passing from electronic hydrogen to muonic hydrogen we observe the variation of 
the relative value of the nuclear structure and polarizability effects, the electron vacuum 
polarization corrections and recoil contributions to the fine and hyperfine structure of the 
energy spectrum llH6j . Muonic atoms represent a unique laboratory for the determination of 
the nuclear properties. The experimental investigation of the (2P — 2S) Lamb shift in light 
muonic atoms (muonic hydrogen, muonic deuterium, muonic helium ions) can give more 
precise values of the nuclear charge radii J7H11|. For more than forty years, a measurement 



of muonic hydrogen Lamb shift has been considered one of the fundamental experiments in 
atomic spectroscopy. Recently, the progress in muon beams and laser technology made such 
an experiment feasible. The first successful measurement of the (fip) Lamb shift transition 
energy (2P[ f = 2 - 2S[ / = 1 ) at PSI (Paul Scherrer Institute) produced the result 49881.88 
(76) GHz (206.2949 (32) meV) [12J. It leads to new value of the proton charge radius 
r p = 0.84184(36) (56) fm, where the first and second uncertainties originate respectively from 
the experimental uncertainty of 0.76 GHz and the uncertainty 0.0049 meV in the Lamb shift 
value which is dominated by the pro ton polarizability term. The new value of proton radius 
improves the CODATA value 13[ by an order of magnitude. Another important project 
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which exists now at PSI in the CREMA (Charge Radius Experiment with Muonic Atoms) 



collaboration proposes to measure several transition frequencies between 2S- and 2P-states 
in muonic helium ions {^\He) + , (/i^Pe) 4 " with 50 ppm precision. As a result new values of 
the charge radii of a helion and a-particle with the accuracy 0.0005 fm will be determined. 
The program of the investigation of the energy levels in light muonic atoms suggests that 
the theoretical calculations of fine and hyperfine structure of states with n — 1, 2 will be 
performed with high accuracy. Note that a discrepancy in the new proton charge radius 
and CODATA value induced both a reanalysis of the earlier obtained contributions to the 
observed transition frequency and a study of the hypothetical muon-proton interaction |14l 
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Theoretical investigations of the Lamb shift (2P — 25), fine and hyperfine structure of 
light muonic atoms was performed many years ago in Refs.[l|, Il8l-l23| on the basis of the 
Dirac equation and nonrelativistic three-dimensional method (see other references in review 
articles |l|, |6|). Their calculation took into account different QED corrections with the accu- 
racy 0.01 meV. Recently an approach of [1] was extended to the case of muonic deuterium 
in [2J, |3] where fine and hyperfine structure was analyzed with high accuracy. Different cor- 
rections to fine and hyperfine structure of muonic hydrogen are calculated on the basis of 
three-dimensional method in quantum electrodynamics in [4j, |24J-|28| . The vacuum polar- 



ization effects of order a 5 were considered in [29l-l3l| . In this work we aim to present new 
independent calculation of the Lamb shift (2P — 2S) in muonic deuterium (/xd) with the 
account of contributions of orders of, a 4 , a 5 and a e on the basis of quasipotential method in 



quantum electrodynamics [264281 |32j. We consider such effects of the electron vacuum po- 
larization, recoil and nuclear structure corrections which are crucial to attain high accuracy. 
With the exception of the nuclear structure and polarizability contribution, we calculate all 
corrections in the intervals (2Pi/ 2 — 2S\/2) and (2P 3 / 2 — 2 Pi/ 2) with a precision 0.0001 meV 
and 0.00001 meV correspondingly. Our purpose consists in a recalculation and improvement 
of the earlier obtained results |l|, |2J and derivation the reliable independent estimate for the 
(2Pi/ 2 — 25*1/2) and (2P 3 /2 — 251/2) Lamb shift, which can be used for the comparison with 
forthcoming experimental data. Modern numerical values of fundamental physical constants 
are taken from Ref. [l3|: the electron mass m e = 0.510998910(13) ■ 10 -3 GeV, the muon mass 
m^ = 0.1056583668(38) GeV, the fine structure constant a' 1 = 137.035999084(51) @, 
the deuteron mass ma = 1.875612793(47) GeV. Numerical values of the proton structure 
corrections are obtained with the 2010 year CODATA value for the deuteron charge radius 



2.1424(21) fm and r d = 2.130 ± 0.003 ± 0.009 fm from [34 



II. EFFECTS OF VACUUM POLARIZATION IN THE ONE-PHOTON 

INTERACTION 

Our approach to the investigation of the Lamb shift (2P — 2S) in muonic deuterium is 
based on the use of quasipotential method in quantum electrodynamics [271, |28|, [35|, where 



the two-particle bound state is described by the Schrodinger equation. In perturbation 
theory the basic contribution to the muon-deuteron interaction operator is determined by 
the Breit Hamiltonian J5|, |36j : 

"^^-"t-^-^ + ^I^ + ^j^- (1) 
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FIG. 1: Effects of one-loop and two-loop vacuum polarization in the one-photon interaction. 



where Hq = p 2 /2fi—Za/r, mi, rri2 are the muon and deuteron masses, \x = vn\ni2i '(1711+1712) . 
The deuteron factor Si — because we used further the common definition of the deuteron 
charge radius r\ 
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The wave functions of 2S- and 2P-states are equal to 
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W = fxZa. (2) 



The ratio of the Bohr radius of muonic deuterium to the Compton wavelength of the 
electron m e /W = 0.7, so, the basic contribution of the electron vacuum polarization (VP) 
to the Lamb shift is of order a(Za) 2 (see Fig. 1(a)). Accounting for the modification of the 
Coulomb potential due to the vacuum polarization in the coordinate representation 
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we present equations for the one-loop VP contributions to shifts of 2S-, 2P-states and the 
Lamb shift (2P- 25): 
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AE VP (2P - 25) = 227.6347 meV, 



(5) 
(6) 



where we round for definiteness the number to four decimal digits. The subscript VP desig- 
nates the contribution of electron vacuum polarization. Experimental error in a determina- 
tion of the particle masses and fine structure constant does not influence on the digits given in 
(6). The muon one-loop vacuum polarization correction of order aiZaf is known in analyt- 
ical form [6]. We included corresponding value AEmvp(2P — 2S) = a 5 /i 3 /307rmf = 0.01968 
meV to the total shift in section V (Eqs.(71)-(72)). This result agrees with that in |2|. Two- 
loop vacuum polarization effects in the one-photon interaction are shown in Fig.l(b,c,d). To 
obtain a contribution of the amplitude in Fig. 1(b) to the interaction operator, it is necessary 
to use the following replacement in the photon propagator: 
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FIG. 2: Effects of three-loop vacuum polarization in the one-photon interaction (a,b) and in third 
order perturbation theory (c). G is the reduced Coulomb Green function (33). 



In the coordinate representation a diagram with two sequential loops gives the following 
particle interaction operator: 
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where the subscript (VP — VP) corresponds to two sequential loops in the Feynman ampli- 
tude (Fig. lb). Averaging (8) over the Coulomb wave functions (2), we find the contribution 
to the Lamb shift of order a 2 (Za) 2 : 
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0.2956 meV 



Higher order a 2 (Za) 4 correction is determined by an amplitude with two sequential 
electron (VP) and muon (MVP) loops. Corresponding potential is given by 
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Its contribution to the shift (2P — 2S) is equal to 

AE VP _ MVP (2P - 2S) = 0.0001 meV 
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The two-loop vacuum polarization graphs (Figs.l(b,c,d)) were first calculated in |39l. 140 . 
The particle interaction potential, corresponding to two- loop amplitudes in Fig.l(c,d) with 
second order polarization operator, takes the form [39]: 
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where the subscript (2 — loop VP) corresponds only to two-loop Feynman amplitudes shown 
in Fig.l(c,d), the spectral function 



f(v) = v{(3-v 2 )(l+v 2 ) 



Lin 



l-v\ n _. n-v\ 3, 1 + v, 1+v . 1 + v, ■ 

+ 2L% 2 + - In In In In v 

1 + v) \l + v) 2 1-v 2 1-v 

(13) 



+ 



ip- 2 >(i + -) + T 



, 1 + v 
In h 



-v(3 — v ) In 2w (3 — v ) In?; 



+ - W (5-3^ 2 )}, 



Li 2 (z) is the Euler dilogarithm. The potential AV 2 C _ X vp( r ) gi yes larger contribution as 
compared with (8) both to the hyperfine structure and Lamb shift (2P — 25). In the case 
of the Lamb shift we find the following contribution: 



AE 2 -ioo P vp(2P ~ 25) = 1.3704 meV 
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Changing in (12) the electron mass to the muon mass one can obtain two loop muon vacuum 
polarization correction. It is known in analytical form from the paper [4l| (we present their 
result with five decimal digits): 
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The numerical values of corrections (9), (14) and the desired accuracy of the calculation show 
that it is important to consider three-loop contributions of the electron vacuum polarization 
(see Fig.2). One part of corrections to the potential from the diagrams of three-loop vacuum 
polarization in the one-photon interaction can be derived by means of equations (8)-(12) 
(sequential loops in Fig.2 (a,b)) [28(. Corresponding contributions to the potential and the 
splitting (2P — 25) are given by 
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where subscripts (VP — VP — VP) and (VP — 2 — loop VP) designate only the Feynman 
amplitudes shown in Fig.2 (a,b) respectively. Six order vacuum polarization contributions 
including Eqs.(18)-(19) were obtained in |2|. The contribution of other diagrams correspond- 



30 for 



ing to the three loop contribution in I7 approximation were calculated in Refs.[29j 
muonic hydrogen. We estimated their contribution for the Lamb shift in (fid) from the re- 
sults given in Eqs.(18) and (23) of Ref. [29J; the result is 0.0021 meV. This gives a total three 
loop contribution 0.0060 meV in one-photon interaction, which is included in Table I. Two- 
loop and three-loop vacuum polarization corrections appearing in second order perturbation 
theory, are calculated in the next sections. Our sum of all three loop VP contributions 
0.0086 meV is very close to the total three loop contribution 0.00842 meV given in Table I 
of [3l| with the account of a rounding. 

Additional one-loop vacuum polarization diagram is presented in Fig.3. In the en ergy 
spectrum it gives the correction of fifth order in a (the Wichmann-Kroll correction) |42j, |43 . 




FIG. 3: The Wichmann-Kroll correction. The wave line shows the Coulomb photon. 



The particle interaction potential can be written in this case in the integral form: 
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Exact form of the spectral function f WK is presented in Refs.p, |42j, |43|]. Numerical inte- 
gration in (20) with the wave functions (2) gives the following contribution to the Lamb 
shift: 

AE WK {2P - IS) = -0.0011 meV. (21) 

This agrees well with other calculations J2j, |3l(. The detailed calculation of all three light- 
by-light graphs is presented in 44|. We included in Table I their estimation using (21) and 
the results from 44 . 



III. RELATIVISTIC CORRECTIONS WITH THE VACUUM POLARIZATION 

EFFECTS 



The electron vacuum polarization effects lead not only to corrections in the Coulomb 
potential (3), but also to the modification of other terms of the Breit Hamiltonian (1). The 



one-loop vacuum polarization corrections in the Breit interaction were obtained in 
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where the superscript B designates the Breit interaction. In first order perturbation theory 
(PT) the potentials AV^p(r) give necessary contributions of order a(Za) 4 to the shift 
(2P-2S): 
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FIG. 4: Effects of one-loop and two-loop vacuum polarization in second order perturbation theory 
(SOPT). The dashed line shows the Coulomb photon. G is the reduced Coulomb Green function 
(34). The potentials AV , AVy P and AVy P are determined respectively by relations (1), (3) and 
(22). 



AE* VP (2P - 2S) = -0.0353 meV, 
AE% VP (2P - 2S) = 0.0011 meV, 
AE^ VP (2P - 2S) = 0.0012 meV, 

AEf VP (2P - 2S) = -0.0023 meV. 



(27) 

(28) 

(29) 

(30) 

The potentials AV^ VP , AV^ VP , AV^ VP take into account the recoil effects over the ratio 
m 1 /m 2 . We have included in Table I the summary correction of order a(Za) A , which is 
determined by equations (27)-(30). The next to leading order correction of order a 2 (Za) 4 
appears in the energy spectrum from two-loop modification of the Breit Hamiltonian. We 
consider in the potential the term of the leading order in m 1 /m 2 (the function f(v) is 
determined by expression (13)): 
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Corresponding (2P — 2S) shift is the following: 
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-0.0002 meV. 
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Other two-loop contributions to the Breit potential are omitted because they give the energy 
corrections which lie outside an accuracy of the calculation in this work. 

In second order perturbation theory (SOPT) we have a number of the electron vacuum 
polarization contributions in orders a 2 (Za) 2 and a(Za) 4 , shown in Fig.4 (b,c): 
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FIG. 5: The three-loop vacuum polarization corrections in second order perturbation theory. G is 
the reduced Coulomb Green function. 



The abbreviation SOPT is used further in Table I, II for the contributions obtained in second 
order PT. The second order perturbation theory corrections in the energy spectrum of 
hydrogen-like system are determined by the reduced Coulomb Green function G (RCGF). 
It has a partial wave expansion [451 ]: 
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(34) 



The radial function g n i(r, r') was presented in [45j in a form of the Sturm expansion in the 
Laguerre polynomials. For a calculation of the Lamb shift (2P — 25) in muonic deuterium it 
is convenient to use the compact representation for the RCGF of 25*- and 2P-states, which 
was obtained in [4j, |46 |: 
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max(xi,x 2 ), C = 0.57721566... is the Euler constant. As 
a result the two-loop vacuum polarization contribution to the first term of (33) can be 
presented originally in the integral form (Fig.4(c)). The subsequent numerical integration 
gives the following results: 
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A B ™?(2P) = -H^- f *Q* f *l)*l x (40) 

x / e - x V+— )dx / e~ x y 1+ ^)dx'g 2P {x,x') = -0.0030 meV, 
Jo Jo 

where the superscript (VP, VP) designates the second order PT contribution when each of 
the perturbation potentials contains VP correction. The results (39), (40) agree with the 
calculation in [3l|. Changing one electron VP potential by the muon VP potential we find 
that corresponding correction to the Lamb shift is very small: 

AE V S ^ T VP (2P - 2S) = 0.0001 meV (41) 

The second term in (33) has the similar structure (see Fig.4(b)). A transformation of 
different matrix elements entering in it is carried out with the use of algebraic relations of 
the form: 
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Omitting further details of the calculation of numerous matrix elements in (42), we present 
here the summary numerical contribution from second term in (33) to the shift (2P — 2S): 

ae sop P t( 2P ~ 25) = 0.0530 meV. (43) 

Other contributions of second order PT (see Fig.4(d,e,f)) have the general structure 
similar to Eqs.(39), (40). They appear after the replacements AVy P — > AV B and AVy P — > 
AVypyp in the basic amplitude shown in Fig.4(c). The estimate of this contribution of 
order a 2 (Za) A to the shift (2P — 2S) can be derived if we take into account in the Breit 
potential the leading order term in the ratio mijva^. Its numerical value is 

AE^'pif' AvB (2P - 2S) = 0.0004 meV. (44) 

The two-loop vacuum polarization contribution is determined also by the amplitude in 
Fig. 4(a). To obtain its numerical value in the energy spectrum we have to use Eqs.(3) 
and (22). In the leading order in the ratio nixjm^ we take again the potential (22), which 
leads to very small correction of order a 2 (Za) 4 : 

AEs£'pt*p(2P - 2S) = -0.00001 meV. (45) 

Three-loop vacuum polarization contributions to the energy spectrum in second order 
perturbation theory are presented in Fig. 5. Respective potentials required for their calcula- 
tion are obtained earlier in relations (3), (8), (12). Considering an accuracy of the calculation 
we can restrict our analysis by a shift of 2S"-level, which can be written in the form: 

AE^ P ^ VP (2S) = -t^f jT p(£R f PWV J™ P (C)d( jH dx(l - |) x (46) 
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FIG. 6: The leading order nuclear structure and vacuum polarization corrections. The thick point 
represents the nuclear vertex operator. 
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In third order perturbation theory (TOPT) the three-loop VP contribution to the Lamb 
shift consists of two terms. One part of it is shown in Fig.2(c). This contribution can 
be calculated by means of (3), (35)-(38) [29|, |3l|. We carry out the coordinate integration 
analytically and the integration over three spectral parameters numerically. The result 



AEtopt' VP ( 2P ~ 25) = 0.0001 meV 



(48) 



is in the agreement with [29|, |31 



IV. NUCLEAR STRUCTURE AND VACUUM POLARIZATION EFFECTS 



An influence of nuclear structure on the muon motion in muonic deuterium is determined 
in the leading order by the root mean square (rms) radius of the deuteron (charge radius). 
We present further all charge radius corrections at two values of r^: r^ = 2.1424(21) fm 
(CODATA 2010) and r d = 2.130(9) fm H(Fig.6(a)): 
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At this point and further the subscript str designates the structure correction. The precise 
value of the deuteron charge radius is needed for the interpretation of new data on transitions 
in muonic deuterium atom. 

There are vacuum polarization corrections connected with the deuteron structure in first 
and second orders of perturbation theory (see diagrams in Fig.6(b,c)). The potential corre- 
sponding to the amplitude in Fig. 6(b) can be written as follows: 
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Its contribution to the 2P — 2S Lamb shift is determined by the formula: 
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FIG. 7: Nuclear structure corrections of order (Za) 5 . The thick point is the deuteron vertex 
operator. 




FIG. 8: Nuclear structure and two-loop vacuum polarization effects in one-photon interaction. The 
thick point is the nuclear vertex operator. 
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The contribution of the same order a(Za) A is specified by the amplitude in the second 
order perturbation theory in Fig. 6(c): 



AE^ SOPT (2P - 25) 



/j, 3 a(Za) 



p(m x 



X- 



367T Jl 

■12 + 23&i - 8bl - Ab\ + Ab\ + 46i(3 - 4&i + 2b\) \nb x 



(52) 



-0.020487 • rj meV 



-0.0940(-0.0929) meV, h 



2m e ^ 
1 + -£. 



Factorizing r\ in expressions (49), (51)-(52) we obtain the finite size correction in the form: 



AE str (2P-2S)+AE^(2P-2S)+AE^ SOPT (2P-2S) = -6.10712-rJ = -28.0309(-27.7074) meV. 

(53) 
The next important correction of order (Za) 5 is described by one-loop exchange diagrams 
(Fig. 7). An investigation of elastic contribution to the Lamb shift and the deuteron polar- 
izability contribution was performed in J22J, |47H50|. Recently new detailed calculation of 
nuclear structure and polarizability corrections which improves previous theoretical results 
is presented in 51J. We have included in Table I the value of the (2P — 2S) shift 1.680(16) 
meV from 15111 . 
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FIG. 9: Nuclear structure and two- loop vacuum polarization effects in second order perturbation 
theory. The thick point is the nuclear vertex operator. G is the reduced Coulomb Green function. 



Two-loop vacuum polarization corrections with an account of nuclear structure are pre- 
sented in Fig.8(a,b,c). The interaction operators constructed by means of Eq.(7) are deter- 
mined by integral formulas: 



AV 



VP-VP/ 
str V 



X 



2Za 3 
27tt 2 

2 



/OO 
p(rj)dr) x 



(54) 



IX 



S(r) - 



in 



r (£2 _ ^ 



(« 



4 e -2m e £r 



av^ vp ( 



AZa 3 
~9^ 



f(v)dv 



o 1 



7r8(r) 



4 e -2m e rir 



mt 



r(l 



VI 



(55) 



The sum of corrections from (54) and (55) to the Lamb shift (2P — 2S) is equal to 

-0.0005(0.0005) meV. (56) 



AEY yP (2P - 2S) 



■'str 



•10.5 -lO -5 -^ 



Two-loop vacuum polarization and nuclear structure corrections of order a 2 (Za) 4 in sec- 
ond order PT shown in Fig.9(a,b,c,d), also can be calculated by means of relations discussed 
in section III. The summary shift is equal to 



AE 



vp,vp 

str,SOPT 



(2P-2S) = -9.5-10 



-5 



-0.0004(0.0004) meV. 



(57) 



There exists also the nuclear structure correction of order a(Za) 5 coming from two- 
photon exchange diagrams with the electron vacuum polarization insertion (see Fig. 10). It 
can be calculated as the elastic contribution of order (Za) 5 J50j. However, there is no need 
to calculate it because in this case we have the same cancelation between elastic two-photon 
correction and deuteron excited states correction as for the contribution of order {Zaf 
5JJ . Indeed, using the notations of Ref. 5l| we can present the muon matrix element Pyp 
for nonrelativistic two-photon exchange with an account of the vacuum polarization in the 
form: 



P^ 



VP 



2a 3 
~3tT 



/OO /■ 



dq 



(47T) 



(2tt) 3 gV + 4ml?) E + £- 



2 

2mi 



3 iq(R-R') 



1 + > 



R 



'\2 



(58) 

where R is the position of the proton with respect to nuclear mass center. Integrating (58) 
over q and expanding the resulting expression over small parameter y/2miE\T{, — R'| we 
obtain: 



32cr 



VP 



/OO 
p(£R 



af - 3a| + 6a ? + 6e" a « - 6 



12af 



(59) 
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-2mi-E|R-R' 



/|2 a ? 



4a| + 12a| - 24a € - 24e~ a « + 24 



48a<? 



a^ = 2m e £|R — R'|. 



It follows from (59) that in the leading order in \J1m\ E\R — ~R!\ elastic correction to atomic 

energy is canceled by the deuteron excited states correction (see more detailed discussion in 

511]). An estimation of second term contribution in the square brackets of (59) to the energy 



spectrum can be derived if we take into account that the integral over £ is determined by 
the region near £ ~ 1. Expanding second term in (59) at small ag we obtain (— 7r/240a^). 
Then performing analytical integration over £ and summing over excited deuteron states we 
obtain the contribution to the Lamb shift: 



6E%(2P-2S) 



2^,3j,2 



m\a 



(0) 



1024m P 



-(<P D \R 2 H D R 2 \<f) D ) - -(folRiHutfRilfo) + 
3 5 



(60) 



+ 2 -{4> D \{R i R j - \8 ij R 2 )H D {R i R j - ^SijR 2 )\<f> D ) 



-0.0001 meV, 



where 4>d is the deuteron wave function. We make all integrations in (60) analytically using 



(61) 



the deuteron wave function in the zero-range approximation [52 

Mr) = \ff-e~ Kr , 
V 2tt r 

where k = 0.0457 GeV is the inverse deuteron size. 

Another term in the Lamb shift of order a(Za) 5 is determined by muon-line radiative 
correction to the nuclear size effect. It was obtained in 53] in a suitable form for subsequent 
numerical estimate: 



AE"i Za)5 (2P- 25) = 1.985 



a(Za) 5 /i 3 



9.62 ■ 10- 



0.0044(0.0044) meV. (62) 



There exists also the correction of order a(Za) 5 with muon vacuum polarization (see dia- 
grams in Fig.(10)). Accounting for the result of its calculation from [6j the total coefficient in 
(62) should be changed: 1.985 —¥ 1.485. However, we can consider together with the muon 
VP and nuclear structure amplitudes in Fig. (10) the contribution of the deuteron excited 
states. Calculating this summary contribution by means of equations similar to (58)- (60) 
(see also 51] ) we observe the cancelation of elastic correction and excited states correction. 
Nuclear structure corrections of order [Zaf can be derived with the use of relativistic 
corrections to nonrelativistic wave functions in matrix element (49) [6|, ]22|, ]54J. We present 
here total contribution to the Lamb shift (2P — 2S) including additional state independent 
correction obtained in |22l . |54j : 

( Za Z n Hr 2 ,\(\nnZ n r\+C-- 



<Za) 



AE ( Jr (2P - 2S) 



12 



{' 



-V- Vd 



(In fiZar) + C 



^+o(r 3 >(-> 



(63) 



jrel 



jrel 



jrel jrel ,2n , ,,2/_4 



)} = -21.28 • 10" 4 • r\ + 0.0029 = -0.0069(-0.0068) meV, 



where the quantities J£ e 3 , F^r are written explicitly in [22|, |54 



We have extracted in the 
square brackets the frequently used quantity (main term) for an estimation of the contri- 
bution to the (2P — 25) Lamb shift in hydrogen atom because other corrections are very 
small (near 1 %) and could be safely omitted. In the case of muonic deuterium they give the 
contribution near 25% of the main term and should be taken into account. Separate energy 
shifts for the 25 and 2P states are given in 22|, ]54J. Numerical estimate is obtained on the 
basis of an exponential parametrization for the charge distribution from J22| . 
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FIG. 10: The nuclear structure and electron vacuum polarization effects in the two-photon exchange 
diagrams. The thick point is the nuclear vertex operator. 





FIG. 11: Radiative corrections with the vacuum polarization effects. 

V. RECOIL CORRECTIONS, MUON SELF-ENERGY AND VACUUM 

POLARIZATION EFFECTS 



An investigation of different order corrections to the Lamb shift (2P — 2S) of electronic 
hydrogen has been performed for many years. Modern analysis of the advances in the 
solution of this problem is presented in a review articles p, |43|, |55|, |56[ • The most part of the 
results was obtained in analytical form, so they can be used directly in muonic deuterium 
atom. In this section we analyze different contributions to the energy spectrum of (fid) up 
to the sixth order in a and derive their numerical estimations in the Lamb shift (2P — 2S). 

There are several recoil corrections of different order in a which give important contribu- 
tions in order to attain the necessary accuracy of the calculation. The recoil correction of 
order (Za) 4 [i 3 /m\ to the Lamb shift appears in the matrix element of the Breit potential 
with functions (2). It is calculated for muonic deuterium in [5l 1371] : 



AE rec (2P - 2S) 



fi 3 (Za) 
12ml 



0.0672 meV. 



(64) 



The recoil correction of fifth order in (Za) is determined by the expression 
fx 3 (Za) 5 a, . 1 8, . , ., 1„ 7 2 



55 



AE (Za) 



m\m<i"Kn a 



-Sio In In fc (n, l)--5 m --a n ^ 

6 Za 6 9 6 vn,2 — vn\ 



2 m 1 2 m 2 > 

■^oio\m 2 In m l in — ) 

ft /i 



(65) 
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no 



where \nko(n, I) is the Bethe logarithm: 

\nk (2S) = 2.811769893120563, 
\nk (2P) = -0.030016708630213, 

«. = _ 2 [h,| + (1+ 1 + ... + 1) + 1 _ i_] Sla + 1(< |' 1)(2| - ; 

Eq.(65) gives the following numerical correction to the Lamb shift: 

AE ( *f(2P - 25) = -0.0266 meV. 
The recoil correction of the sixth order in (Zot) was calculated analytically in [32|, |57j-|6 



(66) 
(67) 
(68) 

(69) 



AEf e f(2P - 2S) = {Za) mi (— - 4 In 2^ = 0.0001 meV. 

8m 2 V 6 / 



(70) 



Omitting explicit form of radiative-recoil corrections of orders a(Za) 5 and (Z 2 a)(Za) A from 
Table 9 [6j, we present their numerical value to the Lamb shift (2P—2S) of muonic deuterium 
atom: 

AE rad _ rec (2P - 2S) = -0.0026 meV. (71) 

The energy contributions obtained in p, |6l|, [62] from radiative corrections to the lepton line, 
the Dirac and Pauli form factors and muon vacuum polarization are given by 



AE 



MVPMSE 



(2S) 



a(Za) jj, 



4 ,,3 



57T m 



4 mi 

- In — ; 

3 fi(Za 



---\nk (2S) + — 



38 

45 



(72) 



a 

71 



9 >/ N 3 2l 10 2 2179\ „ /427 In 2 

-C(3) + -7T 2 In 2 - — vr 2 - -^ ) + AnZa 



2 27 

^E M yp^ MSE (2P) 

a mi /3 



648 ) 

a(Za) A /i 3 



57T m 



n 384 2 

--] Ii U2P)- — 
3 0V ; 6fi 



0.7647 meV, 



(73) 



TT" 



Ti 



3tt // \4 SV ; 2 12 144/ 



-0.0100 meV. 



The diagram in Fig. 11(b) with electron loop polarization insertion in the radiative photon 
gives the contribution to the energy spectrum, which can be expressed in terms of the slope 
of the Dirac form factor F[ and the Pauli form factor F 2 J6(: 



AE. 



rad+VP 



nS) 



/i 3 (Za) 



m\ n 3 



Am(F[(0)5 l0 + F 2 (0) 



C, 



ji 



21+1 



Cji = 5, 



m 



'1-5, 



10 j 



kJV + 1) - /(/ + 1) - f ] mi 



(74) 



(75) 



1(1 + 1) n 

Two- loop contribution to the form factors F[(Q) and F 2 (0) was calculated in [63] (see also 

GIB): 

'a\ 2 rl. ,nn 29 , mi 1 _ 395 

.7T, 



miF[(0) 



iln 2 ^-^ln^i + ic(2) 

9 m P " w ■" " 



108 m e 9 



1296 



(76) 
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TABLE I: Lamb shift (2P]/ 2 — 2S , 1 / 2 ) in muonic deuterium atom. 



/2 



Contribution to the splitting 


AE{2P-2S), meV 


Equation, Reference 


1 


2 


3 


VP contribution of order a(Za) 2 
in one-photon interaction 


227.6347 


(6), [2] 


Two-loop VP contribution of order a 2 (Za) 2 
in one-photon interaction 


1.6660 


(9), (14), [2] 


VP and MVP contribution in 
one-photon interaction 


0.0001 


(11) , [2] 


Three-loop VP contribution in 
one-photon interaction 


0.0060 


(17), (18), [29, 31j 


The Wichmann-Kroll correction 


-0.0011 


(21), [2, Ml 


Light-by-light contribution 


0.0001 


[31] 


Relativistic and VP corrections of order 
a(Za) 4 in first order PT 


-0.0353 


(27)-(30), [5] 


Relativistic and two-loop VP 

corrections of order a 2 {Za) A 

in first order PT 


-0.0002 


(32) 


Two-loop VP contribution of order 
o?{Za) 2 in second order PT 


0.1720 


(39)-(41) , [31] 


Relativistic and one-loop VP 

corrections of order a(Za) 

in second order PT 


0.0530 


(43), [5] 


Relativistic and two-loop VP 

corrections of order a 2 (Za) 4 

in second order PT 


0.0004 


(44)-(45) 


Three-loop VP contribution in 
second order PT of order a 3 (Za) 2 


0.0025 


(46)-(47), [31] 


Three-loop VP contribution in 
third order PT of order o?{Za) 2 


0.0001 


(48), [29, 31], 


Nuclear structure contribution of order (Za) 4 


-27.8749 


(49), [2, 6] 


Nuclear structure and polar izability 
contribution of order {Zaf 


1.6800 


[51] 


Nuclear structure and VP contribution 
in I7 interaction of order a(Za) 4 


-0.0620 


(51) 


Nuclear structure and VP contribution 
in second order PT of order a(Za) 4 


-0.0940 


(52) 


Nuclear structure and two-loop VP 
contribution in I7 interaction of order a 2 {Za) A 


-0.0005 


(56) 


Nuclear structure and two-loop VP contribution 
in second order PT of order oP'iZa) 4 


-0.0004 


(57) 
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Table I (continued). 



1 


2 


3 


Nuclear structure and polarizability 

contribution of order a(Za) 5 

with VP correction 


-0.0001 


(60) 


Nuclear structure contribution of order 
a(Za) 5 with muon-line radiative correction 


0.0044 


(62), [61] 


Nuclear structure contribution 
of order (Za) e 


-0.0069 


(63), [22,54] 


Recoil correction of order (Za) 


0.0672 


(64), [5] 


Recoil correction of order (Za) 5 


-0.0266 


(69), [2,6,55] 


Recoil correction of order (Za) 6 


0.0001 


(70), [6] 


Recoil correction to VP of order 
a(Za) 5 (seagull term) 


0.0002 


[5] 


Radiative-recoil corrections 
of orders a(Za) 5 , (Z 2 a)(Za) 4 


-0.0026 


(71), Table 9 [6] 


Muon self-energy and MVP contribution 


-0.7747 


(72)-(73), [2, 6] 


Muon form factor F[(0), ^(0) contributions 
of order a 2 (Za) 4 


-0.0018 


(78), [4,6,63] 


VP correction to muon self-energy 


-0.0047 


(80), [4,6] 


HVP contribution 


0.0129 


[64,65] 


Total contribution 


202.4139 ± 0.0573 (r d = 2.1424(21) fm) 
202.7375 ± 0.2352 (r d = 2.130(9) fm) 



^(0) 



1 , mi 25 7r 2 m e . m 2 , mi m 2 

- In — + - 4^ In — + 3— % 

3 m e 36 4 mi mf m e m\ 



Then the correction to the Lamb shift is equal to 

AE rad+VP (2P - 25) = -0.0018 meV. 



(77) 



(78) 



To estimate the muon self-energy and electron vacuum polarization contribution in 
Fig. 11 (a), we use the relation obtained in |4J: 



AE 



VP 

MSE 



a 



In 



mi 



d>Tmi{ n(Za) 



c 



< ^ n |A • AV£p|Vn > +2 < ip n \AV£ P GA 



Za 



\lpn > 



(79) 
The sum of all matrix elements which appear in Eq.(79) leads to the following shift (2P—2S): 

AE V M P SE {2P - 2S) = -0.0047 meV. (80) 

The hadron vacuum polarization (HVP) contribution which can be taken into account on 
the basis of the numerical result obtained for muonic hydrogen in 64J, |65j is included in 
Table I. The error of the measurement of the cross section a(e + e~ — > 7r + 7r _ ) was decreased 
to a few per cents. So, we estimate in 5% (±0.0006 meV) corresponding theoretical error of 
HVP correction. 
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VI. FINE STRUCTURE OF THE 2P-STATE 



The leading order (Za) A contribution to fine structure is determined by the operator 

AV fs : 

Za 



AV fs 



Amir 3 



1 + + 2o„ ( 1 + — 

m 2 \ m 2 



fLcr, 



(81) 



AV^ S includes the recoil correction and muon anomalous magnetic moment a M correction. 
Fine structure interval (2P 3 / 2 — 2Pi/ 2 ) for muonic deuterium can be written in the form 



66-68 



AE fs = E(2P 3/2 ) - E(2P 1/2 ) 



_ ^{ZaY 
32ml 

a(Za) 6 fx 3 



3/2J 

2?77i / 7771 

1 + + 2o„ ( 1 + — 

m 2 V 777 2 



(82) 



hm^Zaf mi(Za) 



c, 



+ 



32nmi 



, u(Za) 2 1 

In— — + - 

mi 5 



256 64m 2 

+ a(Za) 4 A VP + a 2 (Za) 4 B VP + A stT (Za)% 2 ■ r\. 



This expression includes the relativistic correction of order (Za) 6 , which can be calculated 
on the basis of the Dirac equation, relativistic recoil effects of order m\[Zaf '/m 2 , correction 
of order a(Za) 6 enhanced by the factor ln(Za) [6], a number of terms of fifth and sixth 
order in a which are determined by effects of the vacuum polarization and nuclear structure. 
Recoil correction (— m\{Za) A I32mfy (the Barker-Glover correction [69]) is also taken into 
account in Eq.(82). This is evident from the expansion of first term in (82) over the mass 
ratio mi/m 2 up to second order terms: mi(Za) A (l — m\jm 2 )j32. The contributions to the 
coefficients Ayp and Byp arise in first and second orders of perturbation theory. Numerical 
values of terms in the expression (82), which are presented in analytical form, are quoted 
in Table II for a definiteness with an accuracy 0.00001 meV. Fine structure interval (82) in 
the energy spectrum of electronic hydrogen is considered for a long time as a basic test of 
quantum electrodynamics [6l l55| . 

Fine structure potential with the leading order vacuum polarization and its contribution 
to the coefficient Ayp are given by (|4|): 



AV&(r) 



AM. r p{s)ds[1 + *2! + 2o „ (1 + i=i )]e -*- 

\2-nm\r 6 J\ m 2 m 2 



AE fs = »MZf \ i + 2m 1 + 2a j i nii 



96irml 



m 2 



m 2 



pttH- 



"(l + 2m e sr)(Lcr 1 ), (83) 
1 + 6^ 



l + 21f0 3 



0.00346 meV. 



(84) 

Higher order corrections a 2 (Zct) 4 contributing to a M are taken into account in this ex- 
pression as well as recoil effects. The same order 0(a(Za) 4 ) contribution can be obtained 
in second order perturbation theory in the form: 



AE V^OPT ~ 1728nm 2 



a(Za) 4 u 3 r , . 2?77i 
' 1 + 2a M + (1 + a„) 



m 2 



p(m 



h (l + 2^) 5 



x 



(85) 



2m e ^ (8m e ^ 



W 



w 



11 



4 1 



2m£ 

W 



In 1 + 



2mJC 

W 



+ 3 



0.00229 meV. 
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Let us consider two-loop vacuum polarization contributions in the one-photon interaction 
shown in Fig.l. They give corrections to fine structure splitting of P-wave levels of order 



a 2 (Za) A . In the coordinate representation, the interaction operator has the form [28l. 135 



AVyp_yp{ 



Za 



1 + 2a u 1 + a. 



Am\ 



+ 



2m\m,2 



[La,) x 



(86) 



X (I?) /"p(0^/"7^^[e 2 (l + 2m e er)e- 2m ^-r / 2 (l + 2m e r / r)e- 2m ^]. 



(e - V 2 ) 

Averaging (86) over the wave functions (2), we obtain the following correction to the interval 

(82): 

fi 3 a 2 {Za) A r, n , 2mi 



x 



AE fs 



P(v)dv 



2887r 2 mf 



*"m4 



1+ 2a M + 



m 2 



[l + a L 



1 



(£2 _ v 2 



-2 ^ w 



1 






„2 ° W I 1 
/ (2m e ri 



+ 1) 3 J 



p(0# x 



0.000003 meV. 



(87) 



V7 i -v V VK 

Two-loop vacuum polarization potential and the correction to fine structure (2P3/2 — 2P1/2) 
are given by 

2Za 3 A+2a h 



Av2-loop VPV 



Lxnj 2-loop VP 



3 ?r 2 r 3 

fi 3 a 2 {Za) 4 



Am 2 



1 + a, 



2m\V(i2 Jo 1 — v 2 



1 f(v)dv — ^^~ 



^r 



1 + 



2m P r 



VT 



OCUrJ, (88) 



A8rc 2 m{ 



9 2 



2m.i , 

l + 2a„ + -1 + a, 

m 2 



/(v)dv ( 6 w 7i-« a + 



1-v 2 1 + 



2m» 



0.00002 meV. 



Wv^ 



(89) 

Two-loop vacuum polarization contributions in second order perturbation theory shown 

in Fig.4(a,d-f) (AV B — > AV^ S ), have the same order a 2 (Za) 4 . For their calculation it 

is necessary to employ the modified Coulomb potential by two-loop vacuum polarization 

271 . l28l ]. The amplitude in Fig.4(e-f) gives the following correction of order a 2 (Za) A to fine 

structure splitting: 



L *- Cj 2-loop VP,SOPT 



1 



1 + 



2m e 



2m e 



wVT 



fi 3 a 2 (ZaY 
?>K 2 mim2 

:+4(l + 



2mi 



f{v)dv 



x 



(90) 



2m„ 



wVT 



ln(l 



2m P 



wVT 



1 — V 2 

= 0.000026 meV. 



Two other contributions from amplitudes in Fig.4(a,d) have the similar integral structure. 
Their numerical values are included in Table II. 

There exists also the correction to fine structure splitting due to nuclear structure. In 
17-interaction it is related with the charge form factor of the deuteron. Fine structure 
potential (81) is obtained for the point deuteron. In the case of the deuteron of finite size 
we can express the contribution of nuclear structure to fine structure splitting in terms of 
the charge radius [68] : 



AE 



f9 = ^(Zaf 2M , 2m! 
64m 2 



-^[1 + 



+ 2^(1 H-^i)] 
m 2 m 2 



-0.00028 meV. 



(91) 



Earlier the calculation of the nuclear structure corrections to the energies of P-levels of order 
(Za) 6 was performed in |22l| . Our numerical result (91) for the fine structure splitting agrees 
with the calculation in |22j . 



TABLE II: Fine structure of 2P-state in muonic deuterium atom. 
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Contribution to fine 
splitting AE? 3 


Numerical value 
in meV 


Equation, 
Reference 


Contribution of order (Za) 4 " 

V?{Za) A (, 2mA 

32m'j V m 2 J 


8.83848 


(82), [2, 6] 


Muon AMM contribution 

Wmf "^ V ~r m 2 ) 


0.01957 


(82), [2,6] 


Contribution of order (Za) e 


0.00031 


(82), [2, 6] 


Contribution of order (Za) & mi/m2 


-0.00001 


(82), [2, 6] 


Contribution of order a(Za) 4 
in first order PT (AVy S p ) 


0.00346 


(84) 


Contribution of order a(Za) 4 
in second order PT 
(AV{? P ■ G ■ AVf s ) 


0.00229 


(85) 


Contributioi 

a(Za) 6 fi 3 
32irm'^ 


1 of order a{Z 

mi 5 


af 


-0.00001 


(82), [6] 


VP Contribution from I7 interaction 
of order a 2 (Za) 4 (AVj s p _ VP ) 


0.000003 


(87) 


VP Contribution from I7 interaction 
of order a 2 (Za) 4 (AV 2 f ° loopVP ) 


0.00002 


(89) 


VP Contribution in second 

order PT of order a 2 (Za) 4 

{AV^p-G-AVJp) 


0.000002 


Fig.4(a), AV B ->• AV fs 


VP Contribution in second 

order PT of order a 2 (Za) A 

(AV<? P _ VP ■ G ■ AVf°) 


-0.000001 


Fig.4(d), AV B -> AW S 


VP Contribution in second 

order PT of order a 2 [Za) At 

{^V£ loovyP -G-AVf s ) 


0.000026 


(90), Fig.4(e-f), AV B -»• AV fs 


Nuclear structure correction 
in I7 interaction 


-0.00028 


(91), [22] 


Summary contribution 


8.86386 





VII. SUMMARY AND CONCLUSION 



In this work, various corrections of orders a 3 , a 4 , a 5 and a 6 are calculated to the Lamb 
shift {2P1/2 — 25*1/2) and fine structure splitting (2P 3 / 2 — 2P1/2) in muonic deuterium atom. 
Contrary to earlier performed investigations of the energy spectra of light muonic atoms in 
[ll, 12, |l8|] , we have used the three-dimensional quasipotential approach for the description 
of two-particle bound state. Our analysis of different contributions to the Lamb shift ac- 
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counts for the terms of two groups. First group contains the specific corrections for muonic 
deuterium, connected with the electron vacuum polarization effects, nuclear structure and 
recoil effects in first and second order perturbation theory. As a rule the contributions of this 
group are obtained in integral form over auxiliary parameters and calculated numerically. 
The necessary order corrections of second group include analytical results known from the 
corresponding calculation in the electronic hydrogen Lamb shift. Recent advances in the 
physics of the energy spectra of simple atoms are presented in the review articles p, |55|, |56[ 
which we use in this study. Numerical values of all corrections are written in Tables I, II, 
which contain also basic references on the earlier performed investigations (other references 
can be found in Ref . [l|, |2|, m) . We compare our intermediate results for different corrections 
with calculations given in |2J. Most part of the results including the Uehling, Kallen-Sabry, 
Wichmann-Kroll corrections, muon Lamb shift contribution, nuclear size and VP corrections 
and recoil terms agree well. Our results for relativistic contributions to the vacuum polariza- 
tion are in agreement with those obtained in [5j. Second order VP correction (39) and (40) 
agrees with the result of [31] just as the three loop VP contribution which is determined in 
Table I by three lines corresponding to one-photon interaction (0.0060 meV), second order 
PT (0.0025 meV) and third order PT (0.0001 meV). Total numerical value 202.4139 meV of 
the Lamb shift (2P — 2S) in muonic deuterium atom from Table I is in good agreement with 
the theoretical result 202.263 meV obtained in [2j. The difference of our result from Ref. [2J 
is connected with the calculation of new contributions of higher order in a and mxjm^-, 
the proton structure and polarizability correction [5l| and slightly different numerical value 
of the charge radius of the deuteron r d used in this work. Two-loop vacuum polarization 
contribution 0.1720 meV of order a 2 (Za) 2 in second order PT is absent in J2|. The value 
of the charge radius r d = 2.139(3) fm is used in [2J. Fine structure splitting (2P 3 / 2 — 2Pi/ 2 ) 
8.86386 meV in Table II agrees also with the result 8.864 meV from [2|. Recently, improved 
analysis of different corrections to the Lamb shift in (fid) is performed in J3(. Total value 
of the Lamb shift (2Pi/ 2 — 2S1/2) for r^ = 2.130 fm according to Table 4 from [3j amounts 
202.9440 meV. This value exceeds our result 202.7375 meV by 0.2065 meV. In our opinion 
the only two essential differences between our Table I and |3( are related with the Zemach 
correction 0.4329 meV and polarizability correction 1.5 meV [3|. It was shown in 5l| that 



the Zemach correction is canceled by the deuteron excited states contribution. As a result 



the nuclear structure and polarizability contribution is equal to 1.680 meV [5jJ which we 
use in our work. 

As has been mentioned above numerical values of corrections are obtained with an ac- 
curacy 0.0001 meV because certain contributions to the Lamb shift (2P — 2S) of order 
a 6 attain the value of tenth part of fieV. The theoretical error caused by uncertainties 
in fundamental parameters (fine structure constant, particle masses) entering the leading 
order contributions is around 10~ 5 meV. The other part of theoretical error is related to 
the QED corrections of higher order. This part can be estimated from the leading contri- 
bution of higher order in a: mia(Za) 6 \n(Za)/7m 3 ps 0.0001 meV. Theoretical uncertainty 
connected with nuclear structure and polarizability contributions is equal to 0.0160 meV 



5JJ. We have also small theoretical uncertainty determined by HVP contribution which we 
estimate in 5% (±0.0006 meV). This estimation is based on the experimental uncertainty 
in the cross section of e + e~ annihilation into hadrons. The rounding errors can amount to 
0.0001^-0.0002 meV. Finally, the biggest theoretical error ±0.0550 meV (for r d = 0.1424(21) 
fm) is related to the uncertainty of the deuteron charge radius. Thereby, the total theoretical 
error of the calculation is equal to ±0.0573 meV. To obtain this estimate we add the above 
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mentioned uncertainties in quadrature. 

Let us summarize the basic particularities of the Lamb shift calculation performed above. 

1. Numerical value of specific parameter m e /fiZa = 0.7 in muonic deuterium atom 
is sufficiently large, so the electron vacuum polarization effects play essential role in the 
interaction of the bound particles. We have considered the one-loop, two-loop and three- 
loop VP contributions to the Lamb shift (2Pi/ 2 — 251/2) • A number of important vacuum 
polarization contributions from 17-interaction agrees with the results obtained in |2|, I29N31 . 



2. Nuclear structure effects are expressed in the Lamb shift of muonic deuterium atom 
in terms of the deuteron charge radius r^. We analyze complex effects due to nuclear 
structure and vacuum polarization in first and second orders of perturbation theory. The 
elastic nuclear structure contribution from two-photon exchange amplitudes is canceled by 
the deuteron polarizability correction J5lj. 

3. Nuclear structure and polarizability effects give the largest theoretical uncertainty in 
the total value of the Lamb shift (2P — 2S). It is useful to express the final theoretical value 
of the (2P - 2S) Lamb shift in the form AE Ls (2P - 25) = (230.4511 - 6.108485 ■ r 2 d ) meV 
with the value of the deuteron charge radius defined in fm. Then, comparing this expression 
with the experimental value of the Lamb shift measured with the precision 0.01 meV (50 
ppm) we can obtain more accurate value of r^ with an accuracy 0.0005 fm. 
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